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Abstract 

The A^-state chiral Potts model in lattice statistical mechanics can be 
obtained as a "descendant" of the six-vertex model, via an intermediate "Q" 
or "r2(t5)" model. Here we generalize this to obtain a column-inhomogeneous 
T2(tq) model, and derive the functional relations satisfied by its row-to-row 
transfer matrix. We do not need the usual chiral Potts relations between 
the A^th powers of the rapidity parameters cip, bp, Cp, dp of each column. This 
enables us to readily consider the case of fixed-spin boundary conditions on the 
left and right-most columns. We thereby re-derive the simple direct product 
structure of the transfer matrix eigenvalues of this model, which is closely 
related to the superintegrable chiral Potts model with fixed-spin boundary 
conditions. 

Introduction 

In a remarkable paper, ^ Bazhanov and Stroganov showed in 1990 how the 
recently-discovered solvable chiral Potts model could be obtained from the 
six- vertex model by a two-stage process. First one looked for a "Q" or T2{tq) 
model whose column-to-column transfer matrix commuted with that of the 
six-vertex model. This turned out to be a spin model on the square lattice, 
each spin taking a given number N of values. Then one looked for a third 
model whose row-to-row transfer matrix commuted with that of the T2{tq) 
model. This was the A^-state chiral Potts model. Some of this working was 
re-presented and extended by Baxter, Bazhanov and Perk. [21 IH] 

Here we focus attention on the square lattice C of L columns. With row 
i we associate a horizontal "rapidity" qi. With column j we associate two 
successive vertical rapidities p2j-i,P2j, as in Figure ^ (except for the initial 
six-vertex model, which has only one rapidity line per column, as we mention 
below) . 

The p and q rapidities are of different types: g is a six-vertex model rapid- 
ity, specified by a single complex variable tq, while p is a chiral Potts model 
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Figure 1: The square lattice £ of L columns with cyclic boundary condi- 
tions, showing the horizontal rapidity lines qi,q2, . . . and the vertical rapidity lines 

Pl,P2, ■ ■■,P2L ■ 

rapidity, specified by four homogeneous variables ap,bp,Cp,dp. Because the 
aim of the earlier papers was to extablish a link between the six-vertex and 
homogeneous chiral Potts models, the chiral Potts conditions 

dp ~\~ k bp = k dp , k dp ~\~ hp ^ k Cp , ( -1 ) 

were immediately introduced. Here k, k' are fixed constants (the same for all 
sites and edges of the lattice), satisfying 

e + k''^ = l . (2) 

The point of this paper is to emphasize that these conditions are not needed 
at the first step of the procedure. Even without them, and with dp,bp,Cp,dp 
allowed to take arbitrary values for each of the 2L vertical rapidity lines of the 
lattice, it is still true that the column-to-column transfer matrix (the "Q" ma- 
trix) of the T2{tq) model commutes with that of the original six- vertex model 
(though not with one another). Further, the row-to-row transfer matrices of 
two models T2{tq), T2{tq) (with different horizontal rapidity variables tq, tq 
but the same vertical p-rapidities) commute. These T2{tq) and Tj{tq) matri- 
ces satisfy straightforward generalizations of the functional relations (4.27) of 
Ref. (21 

We shall also show, still without the conditions (^, that we can define a 
chiral Potts model that is related to the T2{tq) model by appropriate general- 
izations of the transfer matrix functional relations of Ref. [21 • It is, however, 
in general inhomogeneous, its Boltzmann weights being of the usual form, 
but the parameters therein being related in a rather complicated algebraic 
manner both to the p variables and to the tq variable. Further, its row-to-row 
transfer matrices do not in general commute with one another, except in the 
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particular combination 

T{u;'x,Lo^y)f{uj''x,Jy) 

defined in equations (|^^ - ((2HI), - (|^T|) below. ^ (Here i,j,k,l are arbitrary 
integers.) Such a two-row transfer matrix is basically that of the "superinte- 
grable" chiral Potts model. [51 

Finally, we shall remark that this general inhomogeneous model includes as 
a special case the homogeneous T2{tg) model with closed (fixed-spin) boundary 
conditions. The functional relations between the Tj{tq) matrices then simplify, 
and the eigenvalue spectrum is that of a direct product of L single-spin ma- 
trices. This agrees with the properties of such a superintegrable chiral Potts 
model that we observed in [Hj. 

We have written the conditions ((3) down because they are so usually 
associated with the chiral Potts model. ^ Here we never use them. The rapidity 
p = {ap,bp,Cp,dp} has value p{m) = {ap^^),bp^rn),Cp(^rn),dp{m)} on vertical 
rapidity line m, for m = 1,2, ... , 2L. There is no restriction on the complex 
numbers ap(m) , ^p(m) > Cp(m) , <ip(m) • 

The six-vertex model 

We start by defining a six-vertex model in a particular field. For this model 
the doubled vertical rapidity lines pi, ■ ■ ■ ,P2L in Figure ^ should be replaced 
by single "type g" rapidity lines ri, . . . , r^,. 

Associate a spin cxj with each site i of the square lattice C, and allow ai to 
take successive integer values, say 1,2, . . . , N . These can be extended to 
all integer values with the modular N convention = cjj + N. Not all values 
are allowed: vertically adjacent spins aj,ak, with k above j as in Figure ^ 
must satisfy the adjacency rule: 

(Tfc = (Tj or aj — 1 , mod N . (3) 

and horizontally adjacent spins (Ti,aj, with j to the right of i, must satisfy 

aj = Ui or ai + 1 , mod N . (4) 

A typical face k, I, with the corner sites k, I arranged anti-clockwise 
form the bottom- left of the square lattice is shown in Figure ^ With each face 
i,j,k,l associate a Boltzmann weight function WQv{cri,aj,ak,cri). 11 ai = a 
is fixed, then there are only six possible choices of the other three spins, as 
shown in Figure [21 We define the corresponding Boltzmann weights to be 

Wev{a,a,a,a) = cat — 1 , Wev{a,a+l,a,a — l) = ujt — 1 , 
WQv{a, a+1, a+1, a) = t — 1 , WQv{a,a,a — l,a — l) = uj{t — 1) , (5) 
WQv{a, a+l, a, a) = uj — 1 , WQy{a,a,a,a — l) = (w — l)t, 

-"^The T{x,y) of this paper generalizes the T of Ref. W, but with x and y interchanged. 
^They ensure that the column-to-column transfer matrices of the T2{tq) model, i.e. the Q model, 
commute. 
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Figure 2: The six spin configurations of tlie six-vertex model, 
for all integers a. Here 

u; = e^-/^ (6) 

and t is a free parameter. For all other (non-allowed) values of o, b, c, d we 
take WQv{cL,b,c,d) to be zero. We may exhibit the t dependence by writing 
WQv{a, b, c, d) as l^^gy (^l^) b-, c, d). 
The partition function is 

Z = X! n ^6y(crj,cri,o-fc,cr«) ■ (7) 

ijkl 

Here the product is over all faces (i, j, /c, I) of the lattice, with cyclic (toroidal) 
boundary conditions. The outer sum is over all values of all the spins. 

One can regain the usual arrow picture of the six-vertex model by drawing 
arrows on the edges of the dual lattice, pointing to the left or up if the spins on 
either side are equal, to the right or down else. Then the six spin configurations 
of Figure 121 become those of Fig. 8.2 of Ref. 0, where two arrows point into 
each vertex and two point out. Note that the weights © are not those of 
the usual zero-field six- vertex model, since the weights of the third and fourth 
configurations are unequal. A field with weight uj^^'^ has been applied to the 
third and fourth weights.'^ Apart from this field, the \,v of eqn (9.2.3) of [S] 
are related to our present variables uj by e~^^ = lo, = t. 

With each horizontal (vertical) rapidity Qi (rj) we associate a parameter 
tq. ■ Then the model is solvable if for each face 

t = tJU , (8) 

q = Qi being the horizontal rapidity and r = rj the vertical. This is in part 
because the function Wqv satisfies the star-triangle relation 

^6v{tg\b, c, g, a)WQv{tr\a, g, e, DWeviU/tqlg, c, d, e) = 

9 

a, b, g, f)WQvitr\b, c, d, g)W6vitq\g, d, e, /) (9) 

9 

for all tq, tr and all values of the external spins a, b, c, d, e, f. This relation is 
depicted graphically in Figure 01 provided we take Wi , W2 , W3 therein to be 
Wevjtq), Wevitr), WeviU/tq) . 

■^Pasquier and Saleur consider the hamiltonian associated with the six-vertex model in this 
special field. |7| 
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Figure 3: The generalized star-triangle relation 
Transfer matrices 

In this A^-state spin formulation, the row-to-row transfer matrix of the six 
vertex model is an A^'^ by A^^ matrix Uev: with entries 

L 

P(iv]a,b = Y[Weviaj,aj+i,bj+i,bj) , (10) 
i=i 

where a = {oi, . . . , a^} is the set of spins in the lower of two successive rows 
of the lattice C, and b = {bi, . . . , 6^} is the set of spins in the row immediately 
above. It depends on t, so can be written as UQv{t)- 

Similarly, the column-to-column transfer matrix is Vqv = V(ov{t), where 

M 

[^<iv]a,b = YlWQviaj,bj,bj+i,aj+i) , (11) 
i=i 

and a = {ai, . . . , cm} is the set of spins in one column, b = {bi, . . . , 6m}) is 
the set of spins one column to the right, and M is the number of rows of the 
lattice. 

Regarding the spins b,c,e,f as fixed, the star-triangle relation © can be 
viewed as the element (a, d) of an by matrix relation. It involves matrices 
with entries Ws{a,c,d,e), W3{a,b,d, f). Provided these are invertible (which 
they usually are), the relation ensures that [HI §9.6] 

U6v{td)U6v{tr) = U6v{tr)U6v{tg), (12) 

i.e. the row-to-row transfer matrices commute, for all choices of tq,tr- 

Similarly, regarding a, c, d, f as fixed and each side of © as the element 
(b, e) of an A^ by A^ matrix, it also implies that 

V6v{tg)V6v{tr) = VeviU) Vevitg) , (13) 
so the column-to-column transfer matrices also commute with one another. 

The r2{tq) model 

The T2{tq) model is also an A^-state model on the square lattice C, but now 
the spins only need to satisfy the vertical adjacency rule (jH]). The horizontal 
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rule (0J) is not imposed. The partition function is again given by but with 
Wqv replaced by Wr- 

The vertical adjacency rule means that the Boltzmann weight function 
Wria, b, c, d) is zero unless a — d = or 1 (mod N) and 6 — c = or 1 (mod 
A'^). If these constraints are satisfied, then 

Wr{a,b,c,d) = Wr{tq\a,b,c,d) = 

^ LO^^^-'^-Lot^r-^-^Fpgia - d, m)Fp,,{b - c,m) , (14) 

m 

where 

Fpg(0,0) = 1, Fpg(0,l) = -ojcptq/bp, 

Fpg(l,0) = dp/bp, Fpg(l,l) = -uoap/bp . (15) 

This function Wr is the multiplicand of eqn. (3.44a) of Ref. [2]. It is linear in 
the rapidity variable tq and is the Boltzmann weight of two triangles {o, d, m}, 
{5, c, m}, summed over the common spin m, with value or 1, as represented 
in Figure^ The triangles have weights Fpq{a — d,m), Fp'g{b — c,m). There 
are also edge weights w'"'^, u;~"*^, {—iotq)°'~'^, and a site weight (— wtg)"™. 

Here p denotes the four complex variables ap,bp, Cp, dp. Auxiliary variables 
that we shall use are 

Xp — dpi dp , yp — bp/ Cp , tp — XpUp , fip — dp/ Cp . (16) 

Similarly for p'. Throughout this paper we impose no restrictions on ap, bp,Cp, d.^ 
(or ttpi, bpi, Cpi, dpi). They are independent variables. 

In H14|) p and p' are the values of p for the particular face of the lattice 
under consideration. If the face is between spin columns J and J + 1, then 
P = P2J-1 and p' = P2J. 

We see the reason for the doubling of the vertical rapidity lines in Figure 
n The odd rapidities pi,p3, . . . ,P2L-1 are those of triangles such as the one 
on the left in Figure 01 with weight Fpg. The even rapidities P2,P4, ■ ■ ■ ,P2L 
are those of triangles on the right in Figure 01 with weight Fp'g. 



d c 
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Figure 4: The Boltzmann weight Wr{a, b, c, d) of the T2 model as that of two oscu- 
lating triangles. If a, d are in column J and 6, c are in column J + 1, then p = p2j-i 
and p' = P2J. 
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The remarkable feature of the T2(tg) model is that W^- and Wqv together 
satisfy a second star-triangle relation: 

Wr{tg\b, c, g, a)Wr{tr\a, g, e, f)W6v{tr/tq\g, c, d, e) = 

9 

E WGv{tr/tg\a, b, g, f)Writr\b, c, d, g)Wr{tg\g, d, e, /) (17) 
g 

for all tq,tr- 

This can be obtained from Q simply by replacing the first two Wqv 
functions on each side by Wr- It is also represented by Figure 01 but now 
Wi,W2, W3 therein should be replaced by WV(tg), Wr{tr), W^vitr/tq). 

We can define row-to-row and column-to-column transfer matrices T2{tq), 
Vr{p,p') for the T2 model by replacing Wqv in lfTn|) and (fTT|) by Wr- Then 
((T7|) implies that 

T2{tq)T2{tr) = T2{tr)T2{tq) , (18) 

i.e. the row-to-row transfer matrices commute for all tq. 
Also, from (|T7|). 

V6v{t)Vr{p,p') = Vr{p,p')Vev{t), (19) 

for all choices of t,p,p',tq in Q, (|14j) . (|15|) . The column-to-column transfer 
matrices of the six-vertex model therefore commutes with that of any par- 
ticular T2 model. This was the starting-point of Bazhanov and Stroganov's 
derivation^. Note that it does not imply that the column-to-column transfer 
matrices of two different T2 models commute. This is because when uj has the 
particular "root of unity" value © the eigenvalues of the six-vertex model 
are degenerate. 

We shall not consider column-to-column transfer matrices any further 
herein. All the transfer matrices we shall write down in subsequent equations 
will be row-to-row matrices of dimension A^'^ by A^'^. The vertical p-rapidities 
are to be regarded as constants and the horizontal g-rapidity parameters tq 
as variables, in general complex. 

We note in passing that a useful check on both star-triangle relations is 
provided by noting from Q that 

Wevma,b,c,d) = 6{a,c) . 

When tq = tr, it immediately follows that g = d {g = a) on the LHS (RHS) 
of each relation, and that both are trivially satisfied. 

Evaluating (|14j) for the four values of a — d and 6 — c, we find that in each 
case it is linear in the variable tq, so from ()1U|) the matrix T2{tq) is a polynomial 
in tq of degree at most L (usually it is of degree L). From the commutation 
relation ()18() (assuming, as seems to be the case, that the eigenvalues are not 
identically degenerate), the eigenvectors of T2{tq) must be independent of tq. 
The eigenvalues are therefore also polynomials in tq of degree L. 
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The T2, T relation 



The commutation relation H18() is true for any two T2 models with the same 
vertical rapidities pi,p2, ■ ■ ■ ,P2L and different horizontal rapidities tq,tr- Here 
Pj is short-hand for the set {ap^ , bp. , Cp^ , dp^}, so there are actually 8L complex 
numbers specifying the vertical rapidities. We emphasize that there are no 
constraints on these numbers. They can all be chosen independently and H18() 
will still be satisfied. 

The object of this paper is to generalize the transfer matrix functional 
relations of Ref. 2] to this arbitrary inhomogeneous model. We start with the 
T2,T relation of section 4 therein. All equation numbers herein that contain 
a decimal point, e.g. (4.10), are references to equations of Ref. [2j. 

Without loss of generality, we can take A; = in Then (4.4) becomes 

f/-! \ ^ , ,m'd—ma/ , ,j. \a—d—m,' 

d 

Fpq{a-d,m')Fpiq{a-d,m) gj{d) (20) 

and (4.9) is 

E [G.j{a)U,m'{-rj+ir' = g'j{a){-rjr ■ (21) 

m'=0 

Here m, ml take the values 0, 1 and the sum in (|2()j) is over the allowed values 
a, a — 1 of the spin d. 



d 




Figure 5: The sites and faces of C involved in the working of section 4 of If 
are in column J, then p' = P2J-2 and p = P2j-i- 



The RHS of (|20)) is the Boltzmann weight of two successive triangles of the 
T2 model, as shown in Figure |S1 with associated edge and site weights, and an 
additional site weight g,j{d). The spins a, d are in column J ( J = 1, . . . , L) of 
the lattice shown in Figure ^ so the p,p' here are the rapidities P2j-i-,P2,j-2'- 

P = P2J-1 , p' = P2J-2 ■ (22) 
With this identification, (4.11) follows and (4.12) is 

"p2J-2 "P2J-2'J ^ "paj-rg >2J-l'J+l _ -I fr,o\ 



^P2J^2''q '^P2J-2' J >2,/-l ''P2J-1 J+1 
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for J = 1,. . . ,L, taking pq = p2L- This is the condition that the function 
gj{a) be periodic of period N, i.e. gj{a + N) = gj{a). We need this because 
the spins in the T2 model only take A'' values and are always to be interpreted 
as integers modulo N. 

Given , we can solve the bilinear relation (|2,''{|) successively for r^^, . . . , 
rj^, r^_^_-^. Since tl+i = ri, this gives a quadratic relation for , and hence 
for all the . In [2], we then used the fact that we were taking the vertical 
rapidities to be those of the usual chiral Potts model to obtain the two explicit 
solutions (4.13) for the rj. 

Here we can no longer do this: all we can say is that the sequence , . . . 
is one of the two solutions of (|28|) .^ 

We can still carry on with the rest of the working. In place of (4.13) we 

set 

rj = w^-'^^-i xj_i , yj_i = tg/xj^i, (24) 

where is an integer that does not enter the relations H23|). For given 

vertical rapidities pi, ■ ■ ■ ,pl, we take tg,xi, ... ,XL,yi, ■■■ ,yL ^ohe fixed, sat- 
isfying (^5]) and PH) . Then we allow to take any set of integer 
values. 

Equations (4.15) - (4.18) still follow, provided we replace uj~^-'ag/dq by 
Lo'^-'xj, u~^-'Cq/bq by uj~^-'/yj, and similarly with J replaced by J — 1. 
Then in place of (4.19) we obtain 

9j{a) = ypyp'Wj-i{a-l3j-i\ujxj^i,yj^i)Wj{a-l3j\uJXj,yj) , 
9j{a) = — Wj-i{a-l3j-i\xj-i,yj-i)Wj{a-f3j\xj,yj), 

Xpi - Xj-l 



9j{a) = ^^-rZ. Wj-i{a-l3j-i\ujxj-i,ujyj-i) 

X Wj{a-Pj\ujxj,ujyj) . 

Here p = P2J-1, p' = P2J-2 and Xp, yp, tp, Xp',ypi,tp' are defined by ((THl) . The 
functions W, W are given by 



Wj{a\xj, yj) = n '^r'^-^ ^'y-^ly^ , (25) 

= 1 "P2J-1 ^ Cp^j-i^J 



Wjia\xj,yj) = ^'dp^fJ . (26) 

Note the distinction between x and y with a p or p' suffix, and x and y 
with an integer J or J — 1 suffix. The former are defined by ()16l) and are 
vertical rapidity variables. The latter are defined by (^l)) . In fact, xj,yj are 
generalizations of the Xq,yq of [2], so can be thought of as variables", but 
they also depend via (|23|) on all the vertical rapidities. 



*They do have some interesting properties, as in equation l|56|l . 
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Thes functions are generalizations of the chiral Potts model edge-weight 
functions W, W [41 eqns. 2 and 3], ^ eqns. (2.4) and (2.5)]. Their definitions 
can be extended to negative integers a in the usual way: 

a 

= n 

1=1 i=a+l 

for any Sj. However, they do not satisfy the usual periodicity conditions 
W{a + N) = W{a), W{a + N) = W{a). Instead they satisfy the weaker 

condition 

Wjja + N\xj, yj)Wj-,ib + N\x'j_„y'j_,) ^ ^ ^^^^ 

Wj{a\xj,yj)Wj-i{b\x'j_^,y'j_^) 

where x'j = co'^xj, y'j = to^yj and the integers a,b,i,j are arbitrary. In fact 
this is just the condition (^SJ. It ensures that the functions gj{a), g'j{a), g'j{a) 
are all periodic of period A^. 

The "chiral Potts" transfer matrix 

We now generalize the usual definition of the chiral Potts transfer matrix in 
(2.15a) and define a matrix T(x,y) with entries 

L 

T{x,y)a,p = Y[Wj{aj - Pj\xj,yj)Wj-i{aj - Pj-i\xj^i,yj^i) . (28) 
J=i 

Herex = {xi, . . .,xl}, y = {yi, ■ ■ -^yL}, a = {ai, ...,aL} and /3 = . . .,(3l}. 

Because of incrementing any of the spins ai,...,aL by leaves 

T{x,y)afi unchanged. Thus the rows of the matrix T{x,y) have the same 
modulo- spin invariance property as the rows and columns of the T2-model 
tranfer matrix T2{tq). Restricting each of these spins to values, T2{tq) is a 
square A^^ by A^^ matrix; T{x,y) has A^''" rows. 

The columns of T{x,y), labelled by . . . ,(3l, are slightly more subtle. 
Incrementing any /3j by A^ does change the RHS of (|28|) . but only by multi- 
plying it by a factor independent of ai, . . . ,aL. Further, this factor depends 
on xi, . . . ,yL only via their A^th powers. It follows that T{x,y) has at most 
A^^ linearly independent columns, and numerical calculations strongly suggest 
that in general there are indeed A^^ linearly independent columns. Thus al- 
though we may take T{x, y) to have more than A^^ columns, there is a unique 
A^^ by A'^^ matrix Sij{x, y) such that 

Sij{x,y)T{x,y) = T{Jx,Jy) (29) 

for all integers We formally write Sij{x,y) as 

T{uj'x,uj^y)T{x,y)~'^ . 

With these definitions, eqn. (4.20) of becomes 

T2{tq)T{u}x,y) = c{x,tg)T{x,y) + d{y,tg)T{ujx,ujy) , (30) 
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where 

/ . X _ A (j/p2J-l -^Xj){tp^j -tq) 
C[X,lqj — J^J^ ^ , 

j^;^ yp2j-iyp2j \^P2j ^j) 

= n """""^'f"-""'' . (31) 

j=i yp2.j-i yp2j v^P2j-i y^J 

As in (2.42), define an by A^'^ matrix X with entries 

L 

^-,<x' = ll6iaj,a'j + l) . (32) 
J=i 

This is the operator that shifts all spins in a row by one. It commutes with 
T2(t,): 

XT2{tq) = T2(tq)X . (33) 

Replacing xj, yj in (|25|) . (|26|) by l<j~^xj, ujyj is equivalent to replacing the 
index i by i — 1. It follows that 

T{u;-'x,u;y) = p{x,y) X T{x,y) , (34) 

where 

{yp2j-i - - yp2j/yj) 

One other function that we need is given by the obvious generalization of 
(4.23): 

L 

zitq) = \{^ ^^P2J^l^^P2J{^P2J-l-^q)i'tp2J -k)/{yp2.j-iyp2.jf ■ (36) 

J=l 

It is a polynomial in tq, of degree 2L. 
Then one can verify that 

c{x,uitq)d{y,tq)p{u:x,y) = z{ujtq) (37) 

from which it follows that we write 1)30^ as 

T2itq)T{uJX,y) = _f*'^^'',i, . XT{ujx,uj'^y) + d{y,tq)T{uJx,ujy) . (38) 
d{uj ^y,uj Hq) 

The r2,T relation 

We obtained H38|) by considering an A^-dimensional vector g which is a di- 
rect product of L vectors of dimension N, forming T2{tq)g and finding the 
conditions under which this is the sum of two such direct product vectors. 
There are linearly independent vectors g. The matrix with these columns 
is T{x,y). 

We can also form instead the row- vector g^ T2{tq). Corresponding working 
goes through and we are led to define two more Boltzmann weight edge- 
functions 

wM-j^yj) = n r'7r"J-"r ' 

i=i "P2j/yj ^ ''P2.J 
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WMxj.yj) = n """^-^^^'""^""^"-^ (40) 

and a transfer matrix T(x,y) with entries 

L ^ 

T{x,y)i3a = Y[Wj{Pj - aj\xj,yj)Wj-i{Pj-i - aj\xj-i,yj-i) . (41) 
J=i 

Setting 

"(y^^g) = 11 7 7> (42) 

j=i yp2j-iyp2jiyj -^Xp2j-i) 

we obtain the analogue of (4.21): 

f{x,ujy)T2{tg) = ^(^^g) f(x,a;^j/)X + %,tg)f(x,y) . (43) 

The Tj relations 

Here we extend the T2 matrices to the set ri, . . . , tn+i, where tj = Tj{tq) is a 
polynomial in tq of degree {j — 1)L. We shall not give explicit definitions in 
terms of lattice models, as is done in section 3 of Ref. 0, but will use only 
the relation (|38|) . 

We start by defining 

Ar{x,y) = d{y,tq)d{ujy,ujtq) ■ ■ ■ d{uj''~'^y,u''~^tq) , r>0 

= {d{uj-'y,u;-Hq)d{LO-^y,uj~%)---d{uj'^y,uj''tq)}-' , r < 0(44) 

Dr{x,y) = A{x,y)T{ujx,uj'''y)T{ujx,y)~'^ . (45) 

In particular, DQ{x,y) = I is the identity matrix. From (j^UJ, we expect the 
RHS to exist and be unique. 

The relation ()38() can then be written 

T2{tq) = z{tq)XD^i{x,y) + Di{x,y) . (46) 

The spin-shift operator X commutes with T2{tq) and with all the Dr{x,y). It 
is consistent with (|46|) that all the matrices Dr{uj^x, co^y) commute with T2{tq) 
and with one another, for all r, This is what we observe numerically: we 
shall assume that this is so. 
We now define Tj{tq) to be 

i-i 

^i(*g) = Y.^^k)<^k)---''i^^'\)^^ Dk{x,oj^-^y)D,j^k-i{x,i^^y) (47) 

fc=0 

for j > 0. Then 

ro(tg)=0, n{tq) = l 

and T2{tq) is as given in (|46j) . From this definition it follows that 



T2 



{J-hq)Tj{tq) = z{uJ^'hq)XTj^l{tq)+Tj+l{tq), (48a) 
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Tj{ujtq)T2{tq) = z{ujtq)X Tj^i{uJ%) + Tj + i{tq) , (48b) 
TN+l{tq) = z{tg)XTN-liuJtg)+ag + ag, (48c) 

for j = 1, . . . ,N, where 

N-l 

= n diu;'y,u;%), (49) 

i=0 

N-l 

ag = H z{u:Hg)/d{Jy,Jtg) . (50) 

The two sets of relations (|48ap . (I48bl) are equivalent. Note that a^jOg are 
unchanged by the mapping 

x,y,tg x,ujy,ujtg . 

These equations are the generalizations of (4.27) - (4.29). In deriving them 
we have kept x = {xi, . . . ,xl} fixed and incorporated all multiplications by 
powers of to into y = {yi, . . . ,yL} and therefore tg, so if we write the rhs 
of (|^7|l more explicitly as Tj{x,y), then by Tj{u)^tg) in the above equations 
we mean Tj{x,uj^y). However, (|48a|l or ()48b|l can be used to successively 
form T3{tq), Ti{tg), etc. Since z{tg) and T2{tg) are polynomials in tg of degree 
2L L, respectively, from this construction it follows that each Tj{tq) is also 
a polynomial in tg, of degree {j — 1)L. So Tj{tq) is indeed a single- valued 
function of tg, unchanged by replacing x,y by u;x,a;~^y, and by the choice of 
the solution of (|23|) . 

From (|18)) it follows that all the matrices Tj{tq) commute, for all values 
of tq. There is therefore a similarity transformation, independent of tq, that 
simultaneously diagonalizes all the Tj{tg). Then (|48a|) - (|48cj) become scalar 
functional relations for each eigenvalue, which is also a polynomial in tq. These 
relations define the eigenvalue. There are many solutions, corresponding to 
the different eigenvalues. 

Another way of looking at this is to note that if we replace y in (|,38|1 by 
ioy,Lu'^y, . . . ,u!^~^y, we obtain a total of homogeneous linear equations for 
N unknowns T{u>x, y), . . . , T{u!x, u>^^^y). The determinant of these relations 
must vanish, and that is the relation for the function T2{tq) obtained by elim- 
inating T3(tg), . . . ,r7v+i(tg) from (|48a|) and (|48cj) . or equivalently from ()48b|l 
and (|lHc|) . 

We have derived the hierarchy of relations H48a() - ()48c|) from (|38() . We 
could equally well have derived them from (|4Hj) . 

Calculation of aq + olq 

Since each matrix function Tj(tg) is a polynomial in tq, from (|48c() . the same 
must be true of Og + aq, and it must be of degree at most NL. This is by no 
means obvious: it appears from (jlH]) and that Oq and Oq are each quite 
complicated functions of the solution xi, . . . ,xl of ()23() . The object of this 
section is to unravel this little mystery. 
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From (|16l and (|24|) , using the shorthand notation (|22|) , we can write the 
the condition H23|) as 

We can use this relation to ehminate the factors containing yj in H31|l and 
(|i9|) . Using also jSEl) and ((SOI), it follows that 

L ,.N N ( N _ N \(fN_.N\ 

j=i i/p Vp' yyp I 

j=i fp V vV -^j-i) 

where in the multiplicands we again write p,p' for P2J-I1P2.J-21 respectively. 

Since J = 1, . . . , L and xq = xl, is a set of L equations for Xi , . . . , 
We noted above that it has two solutions. Let the other solution be x'j^ , . . . , . 
Then from (EU 

{Xp,—Xj_i){tq—XpXj){tq—yp,x'j_i){yp-x'j ) = (53) 

{Xp, — x' j_i){tg — Xp x'j ){t^ — y^ Xj_i){yp — xj) . 
This equation can be re-written in the "Wronskian" form: 

i^J-l ~ ^'j-l)i'tq — XpXj){yp — x'j ){tp, — tq) = 

(^!/ - -'j){t^ - yP-^-i){^^' - -'j-i){t^ - 1^) ■ (54) 

We can use (|^T|) to eliminate the factors {t^ — XpX^), {t^ — ypiX^_i), 
leaving 



Up yp' K-^p' -^j-iJ 



p yp' 

l^p l^p'K^J ~ ^ J )\^p' ~ ^ J-l)V-p ~ ) 



Vp Vp' ~ ^ J-iJyVp - ^ J ) 



(55) 



Taking the product over J = 1, . . . , L, the factors {xj[ — x'^), {x^_i — x'j_i) 
cancel, giving 

"<7 = [o-q]' ) (56) 

where [oq]' is defined by but with each xj replaced by x'j. Interchanging 
each xj, x'j, it follows at once that Oq = [ZSg] , so + oiq is unchanged by 
replacing the solution xi, . . . ,xl by the alternative solution x'l, . . . , x^. It is 
therefore a single- valued function of tg. 

Now we look at (|^T|) . considered as a recursion relation giving x^_i in 

N 

= fj/gj (57) 



terms of x^ . Set 
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for J = 1, . . . ,L. Then we can choose the normaUzation so that 

- Vp Vp' JJ-1 - (ig - l^p IJ-p'Xp Xp, )fj + [Hp Hp, Xp, - yp, )tg gj 

N N / N N N N\r , / N N.N N N\ /ro\ 

-Vp Vp'Qj-i = (v - ^J'p iJ'p'Xp )fj + [fip i^p'tq - Vp yp> )gj , (58) 

where agam p = P2j-i, p' = V2J-2- This is a hnear relation for {fj-i, gj~i) 
in terms of (/j,gj). 

With these definitions, we find that the multiphcand in the second equation 
((52) is simply gj^i/gj, so 

L 

"9 = n 9J^i/gj = 9q/9l ■ (59) 
J=i 

Define two-by-two matrices 

4.N ,,N 



^V2J t^P2.J 



and set 



fj 
9J 



O = 

Then (|58|) can be written as 

O-l = A2J-2B2J-lij . (61) 

Since Xq = x^ , it follows that = ^Cl, where 

^0 = A a = A2LB1A2B3 ■ ■ ■ A2L~2B2L-i a (62) 

Thus A is the eigenvalue of C/ = A2LB1 ■ ■ ■ A2l~2B2L-i, is the correspond- 
ing eigenvector and, from (jSU, 

oiq = X . (63) 

Since ?7 is a two-by-two matrix, it has two eigenvalues A and A', corre- 
sponding to the two solutions x and x' of the recurrence relations. However, 
we have just shown that interchanging the solutions replaces otq by Oq, so 

aq = X' . (64) 

Since A -|- A' is the trace of the matrix U, it follows that 

Oiq + aq = Trace {A2LB1A2B3 ■ ■ ■ A2L-2B2L-1) ■ (65) 

We are regarding the vertical rapidity parameters Xp^ , yp^ , Hpi 5 • • • ! -^pl iVpl ' 
fjipj^ as constants and tq as a complex variable, so this is an explicit expression 
for aq + clq that makes it clear that it is indeed a polynomial in . Since 
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^2j-2-S2j-i is linear in t^, this polynomial is of degree not greater than L 
(in general it is of degree L). 

From (jnni); and (|Hn]) it is readily seen that 

AT 

AA' = det = z{uj^tq) = aqlXq , 
1=1 

so we could have obained (|64|) without going through the working from equa- 
tion to We have included that working, partly for completeness, but 
also because it is an elegant example of how in solvable models the algebra 
conspires to produce needed results. 

The Tj{tq) relations (|15a|) - (|15c)) . together with (jHOI) and (jHSl), provide a 
closed set of equations that determine the eigenvalues of the Tj{tq) matrices, 
all quantities being polynomials in the complex variable tq. To use them, there 
is no need to solve the eigenvalue equation (|62|) . which is equivalent to the 
recurrence relation H23|) . We could presumably have obtained these relations 
directly by a "fusion" method, generalizing the definition (3.26) - (3.44) of 
Tj{tq), but this is quite technical. We prefer the present approach, based on 
the equation (|38|l. 

We have assumed that the matrices Dr{uj^x,uj^y) commute with one an- 
other and with T2{tq). This assumption agrees with numerical calculations we 
have performed for N = L = 3, but it can probably be removed. We can cer- 
tainly apply a similarity transformation (independent of tq) that diagonalizes 
'^2{tq) (for all tq). Applying this only to the left of (j^H]) . it becomes a set of 
many equations for each eigenvalue of T2{tq). If we focus on just one eigen- 
value and one such equation, then, as we remarked above, we can obtain N 
relations from it by replacing y, tq by w*, w^tg, for i = 0, . . . , — 1. These are 
homogeneous linear relations for the corresponding elements of T(ujx,uj^y), 
so their determinant must vanish. The resulting determinantal relation must 
be equivalent to (|48a|) - (|48c;) . 

The r, T relations. 

From the definitions (|47|) we can also establish that 

aqTj{tg) + Z{tg) ■ ■ ■ z{uJ^'^tq)X^TN~j{ujHq) = Dj{x,U;~^y)Tiy{tq) 

for j = 0, . . . , N . From (|45jl . there must therefore be a matrix Y{x, y), inde- 
pendent of j, such that 

Aj{x,UJ^^y)T{uJX,UJ^^^y)Y{x,y) = aqTj{tq)+z{tq) ■ ■ ■ z{uj^^'^tq)X^TN-j{ujhq) 

(66) 

for j = 0, . . . , A^. 

These equations have the same structure as the fusion hierarchy of relations 
(3.46), except we still have to identify the matrix Y. We have not fully done 
this, but we can note from ^ that 

Wj{a\uJXj,yj)Wj{-a\xj,yj) = 1 , (67) 
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^Wj{a- c\u;xj,yj) Wj{c-b\xj,yj) = 

c 

N{xj-u:-^yj){yp-uxj){t^ -t^) 

{h-t,){x^ -yf){y^ -x^) '^''^ ^''^ 

where p = P2j-i, and 6afi = 1 if a = 6 to modulo N, else 6a^b = 0- The sum 
is over any N consecutive integer values of c: although the W functions indi- 
vidually are not periodic functions, the products in the above two equations 
are indeed periodic functions of a, b, c, of period N. 

Consider the matrix product T{ujx, y)T{x, y), by which we mean the usual 
sum over the intermediate indices, in this case both being the /3 indices in 
(|28|) . (|41|) . Incrementing any /3j by N multiplies the columns of T(ujx,y) by 
certain factors, but divides the rows of T{x, y) by the same factors, so leaves 
their product unchanged. Thus we can naturally take the intermediate sum 
to be over the values 0, . . . , A^ — 1 (or any set of N successive values) of each 
of /3i, 




«J P2J-1 P2J «J+i 



Figure 6: The Boltzmann weight of a face of the transfer matrix product 
T{ujx,y)f{x,y). 

The matrix product is the Boltzmann weight of two successive rows of 
the lattice. A typical face of this double row is shown in Figure |HJ Suppose 
the two external spins a^+i, are equal. Then from (|67l ). the Boltzmann 
weight factors WjWj factors cancel. The sum over the centre spin c then 
gives the RHS of (|^S)) . which vanishes unless aj = bj, mod N. 

If aj = bj, then the same argument applied to the face to the left tells us 
that aj_i = bj^i, and so on. It follows that 

Tiu;x,y)fix,y) = [nitg] + R] /g{x,y) , (69) 

where i? is a matrix with non-zero elements Rab only when ai ^ bi, . . . , 
aiT^bi, and 

a^^^y) - VV^N{x,-u:-^yj){y,-u:xj){t^-t^) ' ^'^^ 

where each p in the multiplicand is p2j-i- 

This is an "inversion identity": it has the same structure as the j = 1 
case of (|66p . with the matrix T{ujx,y) on the left and the first term on the 
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right being proportional to Ti{tq), i.e. to the identity matrix. We conjecture 
(in agreement with numerical calculations) that the right-hand sides of the 
two relations are in fact the same, to within a scalar factor, in which case 
R = z{tq)XTj\f^i{ujtq) /uq and 

Ai{x,uj'^y)Y{x,y) = ag g{x,y) T{x,y) . (71) 

Then (|66|) becomes 

g{x, y) Aj_i (x, y) T{u;x, uj^~'^y) f{x, y) = 

Tj{tq) + Z{tq) ■ ■ ■ z{J-Hq)X^TN-j{i0Hq)/ag . (72) 

This is the generalization of (3.46). 
Consistency 

An interesting consistency check on ((7^ is provided by post-multiplying (|H5|) 
by T{x,Loy), and pre-multiplying (|43j) by T{ujx,y). The left-hand sides are 
then the same, equating the right-hand sides gives 

./ N X T{ujx, uj'^y)f{x, ujy) + d{y, tq) T{ujx, ujy)f{x, toy) = 

^^"^^"^ ■T{uJX,y)f{x,uj^y)X + diy,tq)T{ujx,y)f{x,y) . (73) 



d{ujy,LJtq) 



We can use ((7^ to express each of the TT products as a sum of Tj 
terms. In fact we get only terms proportional to Ti(tg), Ti{ujtq), TN~i{oJtq) 
and TN-i{uj^tq). There are two terms proportional to each of these four fac- 
tors. Since Ti{tq) = I, we can interchange Ti{tq) with Ti{ijjtq) on the RHS. 
Using only the relations 

/\M-2{x,ujy) = aq/{d{y,tq)d{u~'^y,u)'hq)} , 



d{y,tq) = d{y,tq)g{x,y)/g{x,u;y) 
[ with all the 

terms for each r factor cancel, thereby verifying H73|). 



and the commutation of X with all the Tj matrices, we then find that the two 



The T2{tq) model with open boundaries. 

We return to considering the hierarchy of relations H48a() - H48c() for the Tj{tq) 
functions. 

These relations simplify greatly when we impose fixed-spin boundary con- 
ditions on the left and right sides of the lattice. We can do this by taking 

flpi = dpi = . (74) 

Then Fp^q{\, m) = 0, so the weight function of Figure |3 vanishes for the faces 
between column 1 and column 2 unless a = d. This is equivalent to requiring 
that all the spins in column 1 of Figure ^ be equal. The model is unchanged 
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by incrementing every spin by one, so we can in particular require that every 
spin on column 1 be zero. It is evident from Figure ^ that this is the same as 
requiring that all spins on the left and right boundaries be zero. 
((7^ implies that 

fJ-pi = Xp^ = Z{tq) = = (75) 

so the relations H48a() - H48c() simplify to 

T2{tq)T2{uJtq)---T2{uJ^~\) = a . (76) 

From (|65j) . noting that the second row of the matrix Bi is now zero, 

a = [BiA2B^A^---B2L-iA2L]ii • (77) 

The RHS is a polynomial in of degree L and we noted above that T2{tq) 
and its eigenvalues are polynomials in tq of degree L. Further, when either tq 
is large or small, to leading order T2{tq) is diagonal, with entries 

L 

\{{l-u:^J+^-^^+Hq/y2j-iy2j) 
j=i 

in row and column a = {ai, . . . , ai}- Let the zeros of ((77|) be Si , S2 , ■ ■ •, s^. 
Then it follows that all eigenvalues of T2(tg) are of the form 

L 

Aitq) = {u;'^/Y) l[is,-u;^Hq), (78) 
i=i 

where S1S2 • • • sl = uj^^Y, and 

L 

^ = n yp2j-iyp2j ■ (79) 
j=i 

The 7i, • • • , 7l are integers with values in the range 0, . . . , A'"—!. They sat- 
isfy the condition 71 + • • • + 7l = 0, and it seems from low-temperature expan- 
sions that the full set of N^~^ eigenvalues is obtained by allowing 71, . . . , 7l 
to take all such values (distinct to modulo A^). 

It appears that the other values of 71, . . . , 72, that do not satisfy the sum 
rule also correspond to eigenvalues of T2{tq), provided we generalize the model 
to allow the skewed boundary conditions 

aL+i = ai + r 

in every row of the lattice ( so all spins in column 1 are still the same, as are 
all spins in column L + 1, but now those in the two boundary columns no 
longer need be equal). Then 

71 H h 7L = T . 

We can take r to be an integer in the range 0, . . . , — 1. 
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The eigenvalues therefore have the same simple structure as do direct 
products of L matrices, each of size N by A^. For N = 2 this is the structure 
of the eigenvalues of the Ising model. [0] 

For the Ising model this property follows from Kaufman's solution in terms 
of spinor operators i.e. a Clifford algebra. |1H p. 189] Whether there is 
some generalization of such spinor operators to handle the T2{tq) model with 
open boundaries remains a fascinating speculation. |12j 

The results of this section were obtained in Ref. jS]. There we considered 
the superintegable chiral Potts model and rotated it though 90° to obtain a 
model that is in fact the present TM{tq) model. Then we inverted its row-to- 
row transfer matrix, thereby obtaining the present T2(tg) model. We did in 
fact note in section 7 of ^5 that we could allow the modulus k to be different 
for different rows: this corresponds to our here allowing ap,bp,Cp,dp to all 
vary arbitrarily from column to column. 

Summary 

We have shown that the column- inhomogeneous T2(tg) model is solvable for all 
values of the 8L parameters , 6pj , Cp^ , dp^ , . . . , dpj^ , where apj , bpj , Cpj , dpj 
are associated with the Jth vertical dotted line in Figure ^ They do not 
need to satisfy the "chiral Potts" conditions The model then has the 
unusual property that its row-to-row transfer matrices (with different values of 
tq but the same Op^, . . . ,dp^^) commute, while the column-to-column transfer 
matrices do not. 

Our results (jSHl), (ESI), (|^Ha)) - (|lSc)) . (f7^ generahze the relations (4.20), 
(4.21), (4.27a) - (4.27c), (3.46) of Ref. 2 . The last generalization jZl is 
essentially a conjecture, depending as it does on the identification of 
with the J = 1 case of (|66|). However, it has been tested numerically for 
N = L = 3 with arbitrarily chosen values of the parameters and found to be 
true to the 30 digits of precision used. 

One significant difference from the homogeneous model is that the asso- 
ciated chiral Potts model weights p5|) . (|26)) . (|39|) . (|iO|) depend on tq via the 
solution ri, . . . , r^, of H23() . If we change tg then we change ri, . . . , r^, in a non- 
trivial way, so it makes little sense to combine T(x, y) for one value of tq with 
r(x, y) for another value. It appears that our generalized chiral Potts trans- 
fer matrices T and T, with different values of tq, do not satisfy any general 
commutation relations like (2.31) - (2.33) of Ref . [2!. 

In short, we can generalize the T2{tq) model to arbitrary Up-^, . . . , dp^^, 
but the only chiral Potts model we can correspondingly generalize is the "su- 
perintegrable" model with the alternate row-to-row transfer matrices T(x,y), 
T{x, y) defined above. In each double row tq,ri , . . . ,rj^ must be the same for 
T{x,y) and f{x,y). 

The functional relations (|48aj) - (|48cj) define the eigenvalues of the row- 
to-row transfer matrix T2{tq). For fixed-spin conditions on the left and right 
boundaries these can be solved explicitly, giving the simple "direct product" 
result JZHI)- 
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